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Abstract 



We present a novel model to calculate single-electron states in random 
quantum dot superlattices made of wide-gap semiconductors. The source of 
disorder comes from the random arrangement of the quantum dots (configu- 
rational disorder) as well as spatial inhomogeneities of their shape (morpho- 
logical disorder). Both types of disorder break translational symmetry and 
prevent the formation of minibands, as occurs in regimented arrays of quan- 
tum dots. The model correctly describes channel mixing and broadening of 
allowed energy bands due to elastic scattering by disorder. 
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I. INTRODUCTION 



Latest advances in nanotechnology make it possible to growth quantum dot (QD) super- 
lattices In view of the analogy between atoms and QDs, it is expected that strongly 
confined levels overlap when QDs are closely packed. Although this analogy cannot be 
complete since carriers in QDs are influenced by phonons, defects and interface states, the 
resulting collective levels will depend on the particular arrangement of QDs. In this sense, 
QD arrays grown by molecular beam epitaxy can be completely random [|3],fl[], partially regi- 
mented or may be regularly stacked (high regimentation) |5J. Electronic states in highly 
regimented QDs are adequately described within envelope-function approximation with a 
three dimensional Kronig-Penney model and the occurrence of miniband have been estab- 
lished |J. However, the lack of periodicity observed in random QDs superlattices demands 
different approaches. 

In the present work we introduce a two-dimensional effective-mass model to study the 
effects of scattering by disorder on electron transmission through random quantum-dots 
superlattices. To this end, the Ben Daniel-Duke equation is discretized, boundary conditions 
are discussed and scattering solutions are found by means of the transfer-matrix method 
for any arbitrary array of QDs. The model is worked out in a two-dimensional space for 
computational limitations, although it will be clear that generalization to three dimensions 
is rather straightforward. Finally, we present the numerical results for the transmission 
coefficient through random arrays of coupled QDs and the main conclusions of the work. 

II. MODEL 

We consider the Ben Daniel-Duke equation with constant effective mass m* at the T 
valley for the electron envelope function x{Vi z )- 111 order to find numerically the single 
electron states, the whole sample is divided into three different regions, namely left (I) and 
right (III) contacts and the random QDs (II), where scattering by disorder takes place. 
Figure 1 shows a schematic view of the three spatial regions. We then consider a mesh with 
lattice spacings a y and a z in the y and z directions, respectively. Defining t y = — h 2 / {2m*a 2 ) 
and t z = —% 2 /{2m*a 2 z ) 1 we obtain the following discrete effective-mass equation 

^(Xn+l,m "X.n— l,m) ^(Xn,m+1 Xn,m-l) (U n ,rn 2t z 2tjy)^ n)?Tt E\ n ^ m . (1) 

The potential term U n>m in Eq. (|J) is given by the conduction-band edge energy at the point 
(na y ,ma z ). Therefore, disorder enters the equation through this diagonal term. Contacts 
are characterized by flat band conditions, U(n,m) = 0, in the absence of applied electric 
field. Clearly, the effects of the applied field can be easily taken into account within the 
present approach by adding a linearly varying potential of the form Up(n, m) = —eFma z , 
F being the applied field. 

III. TRANSFER MATRIX FORMALISM 

In order to solve the tight binding-like equation ([!]) we use the transfer matrix method 
based on the solutions calculated for each slide of the system along the z direction. For 



2 



the sake of simplicity in the calculation, we define <\> n (n = 0, 1, . . . , N + 1) as a vector 
whose components are 0™ = Xn,m ( m= l>- • • >M). Here M and N + 1 are the number of mesh 
divisions in the y and z directions, respectively. Thus, Eq. (JJ) can be cast in a more compact 
form 



>n-\ 




(2) 



where X and O are the M x M identity and null matrices respectively. The matrix Ai n 
splits into the form Ai n = TZ n + B n . The diagonal elements of the tridiagonal matrix TZ n are 
{7tn)mm — U n , m — 2t z — 2t y while nonvanishing off-diagonal elements equal t y . The matrix 
B n depends on the boundary conditions to be specified later. 

We can obtain the expression of the envelope function amplitudes in the left contact as 
a function of the amplitudes in the right one 




T (N) ( 4>N \ r {N) _ TT ( t z l (El- M, 



(3) 



where T^ N ' is the transfer matrix for the heterostructure. 



IV. SCATTERING SOLUTIONS 

The envelope functions within the contacts will be determined by the boundary con- 
ditions. These boundary conditions are open in the z direction, and periodic on each 
slide, namely in the y direction. Consequently, all elements of TZ n vanish except (H n )iM — 
{Tin) mi = ty The former conditions imply plane wave solutions in the z axis, and the latter 
yield an energy discretization on y. As a consequence, this discretization results in a num- 
ber of transverse channels equal to the number of points in the transverse mesh direction. 
Considering U n ^ m to be constant at the contacts, that is, considering perfect leads so no 
voltage drop occurs within regions I and III, and setting \n,i = Xn,M+i, a particular solution 
of Eq. (|l|) is given by 

X n ,m = ex P h -j^ m J ex P \iha-zn) 

+ 51 rij-j= exp ~jyj rm J ex P (—ikja z n) , (m, n) G I. (4a) 



3=1 \l^3 

at the left contact while at the right contact we have 

M ^ 1 ( 2nj \ 

X l n , m = hi ~7TT exp ( * ~M~ m I exp ( ik 3 a zn) , (m, n) G III. (4b) 



where the normalization constant is choosen to ensure that all the propagating modes carry 
the same current 

Mi = h. sin2 (w) + h sin2 {k ^ ] 
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and 
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E - 2t 
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These expressions remain valid for an applied electric field F provided the electronic mo- 
mentum kj is substituted within region III in the following way 



kj <-> qj 



cos 



1 

2l 



E + eV - 2t 



2wj 
cos — — — 1 
M 



+ 1 



(5) 



Here V = FL, where L is the lenght of region //. 

The matrices f and t appearing in equations fl4aD and ( [Ibf) are the reflection and trans- 
mission matrices and they are responsible for the channel mixing due to scattering events. 
Thus, ?ij represents the probability amplitude for an electron impinging in channel i to be 
reflected into the chanel j. Note that the solution within region II is unknown. Actually, we 
are not interested in this solution since all the physics of the scattering problem is contained 
in the mixing matrices t and f. In particular, we can compute the conductance. From the 
Landauer-Biittiker formalism j7|, the zero temperature two-leads multichannel conductance 
can be calculated using de Fisher-Lee formula M 



2e 2 

G=^Tr(dt), 



(6) 



where Tr stands for the trace of the matrix. 



V. CONFIGURATIONAL AND MORPHOLOGICAL DISORDER 

In order to describe epitaxially random QD superlattices, we consider they are arranged 
on a nonregular lattice (configurational disorder). But, in addition, there exists another 
source of disorder due to spatial inhomogeneities that make the shape to be slightly different 
(morphological disorder). To mimic both types of disorder we assume an array of rectangular 
QDs, randomly displaced from the regular lattice sites, whose size also change randomly 
from dot to dot (see Fig. [I]). To avoid the profusion of free parameters, we consider that 
the energy of the confinement potential AE C provided by the high bandgap semiconductor 
is the same for every QD. This is not a serious shortcoming since spatial inhomogeneities of 
the conduction-band offset and fluctuations of individual QDs shapes yield essentially the 
same results. In addition, our model could easily deal with nonconstant values of AE C if 
further improvements are required. 

We will separate the effects of configurational and morphological disorder for the sake of 
clarity. Configurational disorder means that every QD shifts its position 5r = (5y, 8z) while 
its size d y x d z remains unchanged (see Fig. |]). Here 5y and 5z are random uncorrelated 
variables with zero mean and distributed according to box probability functions of width 
W y and W z , respectively. To simulate the change of the shape (morphological disorder) we 
consider that the QD is enlarged along the Z axis an amount S( while its center stays on a 
regular lattice (see Fig. |2|), where S( is uniformly distributed with zero mean and width W$. 
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VI. RESULTS 



As a working example, we have performed several numerical calculations in order to 
study the effect of both configurational and morphological disorder over the conductance of 
random QDs superlattices made of GaAs-In-Ga^-As heterojunctions. We have considered 
regimented and disordered 4x4 arrays of QDs to elucidate the effects of randomness. The 
conduction-band offset, AE C , is taken to be 70% of the difference of the gaps AE g in strained 
GaAs-Inj.Gai_.jAs heterojunctions, where AE g = 1.45a; eV. For defmiteness we set x = 0.35 
and consequently AE C = 0.35 eV. In addition, since we are mainly interested in the effects 
of the coupling between the QDs through the high bandgap semiconductor rather in the 
confined levels of individual QD, we have taken m* = 0.067 in units of the bare electron 
mass, corresponding to the embedding semiconductor. Let us mention that the the model 
can be easily generalized to include a different effective mass inside the QDs. The size of 
the regimented QDs is d y x d z with d y = 8.0 nm and d z = 1.6 nm. The separation between 
centers in the regimented array is 14.0 and 6.8 nm along the Y and Z axes, respectively. 
The number of mesh points along the two spatial directions are M = 50 and N + 1 = 39. 

As typical results of our simulations, Fig. |3| shows the conductance versus Fermi energy, 
measured from the conduction-band edge in In x Gai_ x As, in the absence of applied field. 
Solid line corresponds to regimented QDs. The coupling between the QDs splits the energy 
levels and results in the formation of minibands ||. From Fig. |3| we observe the occurrence 
of two well-defined minibands below the barrier when disorder is absent, in agreement with 
previous results ||. Each band is characterized by four main conductance peaks and each 
peak is the convolution of four peaks that cannot be resolved except for the lower one in the 
higher miniband (see inset of Fig. ^). 

Transport though the miniband changes as soon as some degree of randomness is consid- 
ered in the model, as expected. Mean size and separation between QDs as well as fluctuations 
around those values strongly depend on the growth conditions (e.g. growth temperature) 
and subsequent thermal treatments Q. As an example, Fig. |3| shows the conductance for 
configurational (W y = 2.0 nm and W z = 1.2 nm) and morphological (W^ = 0.8 nm) disorder, 
obtained by averaging over 100 realizations of the disorder. As a main point, we notice that 
electronic states in the array of random QDs behaves like an amorphous material since the 
conductance strongly decreases while the allowed energies broadens due to the fluctuations 
of QD energy levels for each realization of the disorder. 

Concerning the effects of a uniform applied electric field, we have computed the con- 
ductance for a given Fermi energy (E = 0.27 eV) as a function of the applied bias V. The 
potential was assumed to drop uniformly across region II. Figure f| shows the conductance 
versus applied bias for regimented as well as random QD superlattices. Regimented QD 
superlattices present three well-defined negative dG/dV regions due to resonant tunneling 
through the QDs. The observed peak-to-valley ratios become worse as soon as disorder is 
included in the calculation, especially in the case of configurational disorder. 

VII. CONCLUSIONS 

In summary, we have presented a method to study electron transport through random 
QD superlattices. The method is based on the transfer matrix formalism applied to the 
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discrete Ben Daniel-Duke Hamiltonian for the electron envelope function. A careful analysis 
of the scattering solutions under appropriate boundary conditions (periodic and open along 
the lateral and longitudinal directions, respectively) allows us to obtain the two-channel 
conductance. For regimented QDs (regular array of QDs) the conductance shows clear sig- 
natures of the miniband structure, as previously predicted by a Kronig-Penney model for 
strongly coupled QDs ||. However, the novelty of the model lies in the fact that random 
QD superlattices can also be studied. To this end, two models of disorder (configurational 
and morphological) have been introduced. Disorder reduces the conductance due to An- 
derson localization of the envelope functions and broadens the allowed energy bands. The 
characteristic G-V presents several regions of negative dG/dV, although the peak-to- valley 
ratios strongly decrease due to disorder. 
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FIGURES 



FIG. 1. Schematic view of the sample. Regions I and III are the electrical leads of the sample 
(contacts) and electrons undergo scattering processes only at region II. 

FIG. 2. Schematic view of configurational and morphological disorder. In the former case the 
QD shift its position an amount 5r while its nominal size dy x d% remains unchanged. In the later 
case, the center of the QD does not change while its size along the Z axis increases or decreases 
an amount <5£. 

FIG. 3. Conductance versus energy for a 2D ordered array of 4 x 4 QDs made of In^Gai-^As 
in GaAs (solid line) with no applied electric field. The inset shows an enlarged view of the lower 
conductance peak of the second miniband. Results are compared to random 4x4 arrays with 
configurational (dotted line) and morphological (dashed line) disorder. 

FIG. 4. Conductance versus applied bias for a 2D ordered array of 4 x 4 QDs made of 
In x Gai_ x As in GaAs (solid line), when the incoming electron energy is E = 0.27 eV. Results 
are compared to random 4x4 arrays with configurational (dotted line) and morphological (dashed 
line) disorder. 
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Fig.l, I. Gomez et al. 
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Fig.3, 1. Gomez et al. 
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Fig.4, 1. Gomez et al. 
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